Flavor Symmetry and Competing Orders in Bilayer Graphene 
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We analyze competition between different ordered states in bilayer graphene (BLG). Combining 
arguments based on SU(4) spin-valley flavor symmetry with a mean field analysis, we identify the 
lowest energy state with the anomalous Hall insulator (AHI). This state is an SU(4) singlet excitonic 
insulator with broken time reversal symmetry, exhibiting quantized Hall effect in the absence of 
external magnetic field. Applied electric field drives an Ising-type phase transition, restoring time 
reversal symmetry. Applied magnetic field drives a transition from the AHI state to a quantum 
Hall ferromagnet state. We estimate energies of these states, taking full account of screening, and 
predict the phase diagram. 

PACS numbers: 



Bilayer graphene (BLG), due to its unique electronic 
structure of a two dimensional quadratically dispersing 
semimetal [H, offers an entirely new setting for investi- 
gating many-body phenomena. The density of states in 
BLG does not vanish at charge neutrality, and thus even 
arbitrarily weak electron interactions can trigger phase 
transitions in charge neutral system. Theory predicts in- 
stabilities to numerous strongly correlated gapped exci- 
tonic states with different spin and valley structure 0-0] , 
as well as gapless nematic states [f|. 

Experiments indicate that the gapped state observed 
in BLG at charge neutrality in quantizing magnetic fields 
@ persists down to low fields, crossing over to another 
gapped state at zero field Q ■ Gapped states in BLG can 
feature new interesting properties. Also, they can lead to 
new applications in nanoelectronics, in particular those 
based on switching between gapped and ungapped states 
in external fields. Hence, clarifying the relation between 
different gapped states and understanding the phase dia- 
gram in external fields is an interesting and timely task. 

Excitonic instability results from the interaction mix- 
ing conduction and valence band states and sponta- 
neously opening a gap between them. A variety of or- 
dered states with different symmetries (ferromagnetic, 
fcrrimagnctic, antifcrromagnetic, ferroelectric, etc) have 
been analyzed using the short-range interaction model 
[3, IH and dynamically screened Coulomb interaction Q . 
Surprisingly, all predicted gapped states appear to be 
degenerate on a mean field level, suggesting an under- 
lying symmetry relation between them. Understanding 
the competition between these states calls for a unifying 
symmetry-based approach which will be outlined below. 

Our analysis indicates that the state realized at charge 
neutrality breaks Z% time reversal symmetry but is in- 
variant under a continuous symmetry group SU(4). This 
state, exhibiting quantized Hall effect in the absence of 
magnetic field |8| , is identified with one of the topological 
insulator states discussed recently, the anomalous Hall in- 
sulator (AHI) [9|. Applied electric field drives an Ising 
type transition from the AHI state to a layer polarized 
gapped state in which time reversal symmetry is restored. 
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FIG. 1: a) Schematic phase diagram for BLG at charge neu- 
trality. Electric and magnetic fields drive transitions from the 
anomalous Hall insulator (AHI) state, realized at low E and 
B, to the layer polarized state (high E) and quantum Hall fer- 
romagnet state (high B). As discussed in the text, ordering in 
these states is described by a 4 x 4 matrix Q = 1 (AHI state), 
Q = f)3 (E field insulator), and Q = 03 or Q — 771.2 (spin or 
valley polarized QHFM state), b) Landau level spectrum of 
the AHI and QHFM states (E — 0) . Note an anomalous Lan- 
dau level in the AHI state that has no particle-hole-symmetric 
counterpart. Particle-hole symmetry is restored in the QHFM 
state, allowing the chemical potential to relax to zero. 

Applied magnetic field drives a transition to a quantum 
Hall ferromagnet state [see Figflja)]. 

As we shall see, these states are governed by SU(4) 
symmetry defined in terms of four spin and valley fla- 
vors. Such symmetry is peculiar to BLG: although the 
2x2 Hamiltonians describing particles in valleys K and 
K' have different chirality, they can be mapped on each 
other by permuting the wavefunction components on sub- 
lattices A and B. Importantly, since the main interac- 
tions are SU(4) invariant, all gapped states can be classi- 
fied according to different subgroups of the SU(4) group 
which describe different patterns of spontaneous symme- 
try breaking. Weaker interactions can then be incorpo- 
rated as external fields in a sigma model framework. As 
we show below, this approach can be successfully im- 
plemented to construct a full phase diagram of BLG in 
external electric and magnetic fields. 

We start with recalling some basic facts about elec- 
tronic structure of BLG. The low-energy electron states 
can be described by a two-component wavefunction tak- 
ing values on the A and B sublattice of the upper and 
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lower layer respectively To analyze the structure of 
the Hamiltonian, it will be convenient to combine the 
spin and valley components in a single eight-component 
wavcfunction 4'a,s,v{ x -)j where a is the sublattice (layer) 
index. We shall use the Pauli matrices in sublattice, spin 
and valley space, denoted below by t%, and rji, respec- 
tively. The low energy non-interacting Hamiltonian may 
then be written as 



H = 



{Px+iPymY 



{px - ipymf 



2m 2m (1) 

where t± = t\ ± ir 2 . Here m = 0.05m e is the effective 
mass, and the parameter e describes the effect of trans- 
verse electric field (see discussion below). 

Because of the presence of 773 in Eq.(fT]), the Hamil- 
tonian is not invariant under rotations of valley compo- 
nents. To bring it to an SU(4) invariant form, we perform 
a unitary transformation on all operators 

1 + m 1 



o = uou\ u 



V3 



-Tl 



(2) 
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This transformation mixes the layer and valley indices of 
the wavefunction i)) a ^ v (x) by interchanging layers in one 
of the valleys. Defining p± = p x ± ip y , the transformed 
non-interacting Hamiltonian takes the compact form 



fr - p + + 

2m 



p 2 ^ 



2m 



(3) 



At e = 0, this Hamiltonian is manifestly invariant under 
SU(4) rotations in the spin/valley flavor space. 

Electron interactions can be described by a many-body 
Hamiltonian written in terms of p q = -0 P V-'p+q (the 
density summed over layers) and A q = ^ p ^p^Vs^p+q 
(the density difference between layers). The interacting 
Hamiltonian, which incorporates a difference between in- 
terlayer and intralayer interaction [10| , can be written as 



H 



2^L^+(9W-q + ^-A q A- q , (4) 
p q 

where V+(q) = 27re 2 / nq is the Coulomb interaction, and 
V- = Tre 2 d/n accounts for the layer polarization energy 
(here d = 3.5A is the BLG layer separation). The pp 
term, which is isotropic in flavor space and thus is SU(4) 
invariant, dominates because d is small compared to 

10k A, 



ao = h n/me 



(5) 



the characteristic lengthscale set by interactions [3|. 
Working in the long wavelength limit, qd -C 1, we shall 
make an approximation V_ •CV+, by suppressing the AA 
term which breaks the SU(4) symmetry, and restoring it 
at a later stage in our analysis. 

Excitonic instability, analyzed in Refs.0-3], results in 
correlated states described by the BCS-like hamiltonian 



P+f- 



P-f^ 



2m 



Ar 3 Q ) tl)p+-V + (p)p p p- p , 
(6) 



with A the gap order parameter and Q a 4 x 4 hermitian 
matrix in flavor space satisfying Q 2 = 1. Mean field 
theory fixes the value of A up to a sign, but leaves 
Q undetermined, since all choices of Q give the same 
energy. 

Since hermitian matrices satisfying Q 2 — 1 have eigen- 
values ±1, all excitonic states can be classified as (a, /3), 
where a and j3 are the numbers of +1 and —1 eigenvalues 
of Q respectively. This gives three types of states: (2, 2), 
(3,1), and (4,0). There is an additional Z 2 symmetry 
associated with the overall sign of Q which is absorbed 
into the sign of A. 

Symmetry of these states depends on the ordering 
type: U(l) x SU(2) for Q in the (2,2) manifold, and 
SU(3) for Q in the (3,1) manifold. The case (4,0) is 
special: Q is a unit matrix, and the SU(4) symmetry 
of the Hamiltonian is not broken. This SU(4) singlet 
state, which polarizes layers by valley, is the valley anti- 
ferromagnet of Ref.0]. Other states discussed in Refs.@- 
4] correspond to continuous symmetry breaking of type 
(2, 2) or (3,1). 

While within the mean field theory all choices of Q 
yield the same energy, this degeneracy can be lifted by 
the effects due to order parameter fluctuations. These 
fluctuations, which tend to suppress ordering, will be 
stronger for states (2,2) and (3,1) which break contin- 
uous symmetry and thus host gapless Goldstone modes 
than for the SU(4) singlet state Q = 1. Hence we expect 
A( 2j 2)j A(3,i) < A(4 ;0 ). Meanwhile, the nematic state 
proposed in 0] has a higher free energy even at mean 
field level (see supplement). We therefore predict that 
BLG should spontaneously order into the Q = 1 state. 

While the Q — 1 state does not break continuous sym- 
metry it nonetheless has interesting properties. In par- 
ticular, the mass term Af 3 in Eq.© breaks the time 
reversal symmetry implemented for the Hamiltonian, 
Eq.(|4]), by T = Kf]\ ®f\® ia 2 , where K represents com- 
plex conjugation. The Q = 1 state is thus a state of the 
type proposed in Q as a condensed matter realization 
of the parity anomaly, called 'anomalous Hall insulator' 
(AHI) in the papers on topological insulators [9(. Spon- 
taneous breaking of time reversal results in a quantized 
Hall effect (QHE) even at zero magnetic field. 

The anomalous Hall effect can be understood in terms 
of edge states. Since the choice of sign for A breaks a Z 2 
symmetry, the AHI state falls into the universality class 
of the Ising model. We therefore expect spontaneous 
formation of domains, with A(x) changing sign across 
the domain boundaries. As shown in [111 ], such domain 
boundaries host topologically protected edge modes. In 



contrast to [11] ]. however, all edge modes have equal chi- 
rality, resulting in charge QHE rather than valley QHE. 

Because domains with opposite signs of A will have op- 
posite o~ xy , the Hall conductance of a macroscopic sample 
will average to a value near zero. At the same time, per- 
colating edge states will contribute to the longitudinal 
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conductance cr xx . 

We now consider applying a transverse electric field to 
BLG, e 7^ in Eq.Q, which imbalances the potentials on 
the two layers. This reduces the symmetry group of the 
Hamiltonian to SU(2) sp i n <£> U(l) va iiey) with generators <7j 
and 773 respectively. We note that electric field favors po- 
larizing the layers by charge, i.e. Q = r) 3 . This state 
has the same symmetries as the Hamiltonian (at e 7^ 0), 
and in particular preserves time reversal symmetry. Thus 
an electric field induces a phase transition in which time 
reversal symmetry is restored, and the anomalous Hall 
conductance disappears. This transition is analogous to 
restoration of Z2 symmetry in a quantum Ising ferromag- 
net upon application of an external magnetic field [121 ]. 

We estimate the critical electric field E c by equating 
e s» A 4-0 — A 2i 2- Here e — ^eEd/n^, where E is the 
electric field that would exist in the absence of graphene, 
and K- describes intrinsic screening. The value of k_ 
can be found by following [l3| and cutting the infrared 
logarithmic divergence by the gap A. Using the value Q 



1.5meV« 



(7) 



gives k_ ~ 3. For an order of magnitude estimate, we use 
A^o instead of A^o — ^2,2, finding E c ~ 2.6K~ 1 meV/A. 

Including the interaction V- , which disfavors layer po- 
larization by a 'capacitor energy', Eq. (J4j , does not change 
the above conclusions. While penalizing the Q = 773 
state may shift the critical value E c , it does not alter 
the qualitative picture discussed above, whereby a time 
reversal symmetry that is broken at E < E c is restored 
at E > E c . Furthermore, the quantitative effect on E c 
should be parametrically small in d/a$. 

Including trigonal warping [l4| in our Hamiltonian sim- 
ilarly does not alter our conclusions. The trigonal warp- 
ing term is not SU(4) invariant. However, it is time rever- 
sal invariant. Thus, when trigonal warping is sufficiently 
weak, it does not alter ordering in the AHI state. 

The flavor symmetry framework introduced above can 
be applied to BLG in a magnetic field. The main effect of 
magnetic field is orbital coupling, described by p = p 
eA in Eq.©, which preserves the SU(4) symmetry of the 
Hamiltonian. This orbital coupling causes the spectrum 
to split into Landau levels [l4| with an energy spacing 
of order huj c , where oj c = eB/mc. Magnetic field also 
couples to spin via Zeeman interaction, characterized by 
Ez = 2[1bB. However, since this interaction is weak, 
Ez <C hu c , we proceed by neglecting it at first. 

We expect the Q = 1 ground state introduced above 
to be robust to the application of a weak magnetic field. 
However, the T non-invariance of the mass term AT3 
means that the Landau level spectrum is not invariant 
under B — > —B [l5j |. In particular, the zeroth Landau 
level forms at energy +A only, and has no counterpart at 
—A. If BLG is to remain at charge neutrality in a mag- 
netic field, the chemical potential p must be pulled up to 
the zeroth Landau level [see Fig[TJb)]. In a system with 



domains, meanwhile, the chemical potential will remain 
fixed, resulting in band bending near domain boundaries. 

At a higher magnetic field the system will transition 
to a quantum Hall ferromagnet (QHFM) state [3], in 
which spontaneous flavor polarization of the zeroth Lan- 
dau level is driven by exchange energy. Within our SU(4) 
symmetry framework, QHFM ordering corresponds to Q 
of type (2,2). A Q 2 ,2 mass term in Eq.((6|) describes gap 
opening in the zeroth Landau level, with spin or valley 
polarization controlled by the choice of Q. We now re- 
view some basic properties of the QHFM state, taking 
full account of screening, which was neglected by 16 1. 

When evaluating the energy gain from QHFM order- 
ing, it is essential to take into account screening aris- 
ing from Landau level polarization. This is so because 
the characteristic energy scale for interactions within the 
zeroth Landau level e 2 /kIb greatly exceeds the Landau 
level spacing ftw c for all experimentally relevant magnetic 
fields (here I b = y/hc/eB is the magnetic length) . As we 
shall see, the QHFM energy is controlled by the behav- 
ior of the screened Coulomb interaction V+ (q) in the l ong 
wavelength limit qls J$ 1. Applying the result of Ref. 17 1 
in this limit, we find 



V+(ql B < 1) = 



27re 2 /K 
q + \Zq 2 l B 



Z = N^»1, 
a 



(8) 



with N = 4 and A = 0.88 a numerical prefactor. In a wide 
range, Z^ 1 < qls < 1, this expression behaves as l/q 2 , 
giving rise to a logarithmic spatial dependence V+ (r) de- 
scribing field-induced 'unscreening' of the screened po- 
tential discussed in Ref. [l8|. Such unscreening will be 
manifested by any state with a dipole active gap. 

Additional insight into the behavior of the screened 
potential Eq.(j8]) can be obtained by comparing to the 
dynamical polarization function found at B = and in 
the absence of any gap [3J]. Evaluating it at frequency 
u> = u) c , in the limit qls < 1 we find the interaction 
identical to Eq.®, with A = 1. 

Having obtained the screened potential V+ , we proceed 
to calculate the QHFM exchange energy gain per electron 
J(B) in a 'jellium model' [l9| as 



d 2 p 



(l-g(p))V + (p). 



(9) 



Here, N — 4 reflects the number of occupied states in the 
zeroth Landau level, and g(p) is the Fourier transform of 
the exchange hole g(r) = (po(r)po(O)) / (po{0)) 2 , where 
pa is the electron density projected onto the zeroth Lan- 
dau level, and the average is taken with the many body 
wavefunction |*) of the system ((A) = 

When calculating g(r), it is important to note that the 
zeroth Landau level is eightfold degenerate, exhibiting an 
accidental orbital degeneracy in addition to the fourfold 
flavor degeneracy. However, as was shown in [lij . flavor 
polarization is favored over orbital polarization. Since 
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the QHFM exchange energy is the same for any flavor 
polarized state, without loss of generality consider a state 
|\&) polarized in spin. The exchange hole g(r) for such 
a state may be calculated by working in Landau gauge, 
and remembering that single particle states with different 
flavor are orthogonal, whereas single particle states with 
different orbital index or guiding center coordinate are 
not orthogonal. After some algebra, we obtain 



;M = l-e-« 2 / 2 (£ 2 -4) 2 /32, i = r /l B . 



(10) 



We note from Eq.© and Eq.([TU]) that the exchange en- 
ergy comes predominantly from wavevectors qls J$ 1, as 
expected. Using the 1/q 2 potential, we obtain, with log- 
arithmic accuracy, 



J(B) 



2tt e 2 a m(l + Z) 



n\l 2 B 



— hu! r In — - 
8A B 



(11) 



where <& = hc/e is the flux quantum, and B is given 
by the condition Z = 1. This is the natural result for a 
particle interacting with ~ B other particles a distance 
~ B~ l l 2 apart in a logarithmic potential, and differs 
strongly from the E ~ \J~B~ prediction from [16] ■ The dif- 
ference arises because we have taken account of screening. 

The critical magnetic field B c for the field induced 
transition from the AHI state to the QHFM state may 
be estimated by equating J(B) — A^o, which gives 
B c ~ 0.5 T. This is in agreement with recent observa- 
tion of a resistive state in BLG at v = 0, appearing at 
B > 0.1T The resistance of this state was observed 
to increase exponentially with magnetic field, where the 
exponent scaled linearly with Bjk^T. This is consistent 
with our analysis above, assuming that transport in the 
QHFM state is controlled by thermal activation across 
the gap J{B). 

The effect of external electric field on the QHFM state 
can be conveniently analyzed in a 'sigma-model' frame- 
work, where the electric field, the 'capacitance energy' 
V_, and the Zeeman term are treated as perturbations. 
We suppress the unimportant orbital index, and char- 
acterize the flavor structure of the QHFM state by a Q 
matrix of type (2,2). We can write the free energy per 
electron F = j d 2 rf[Q(r)}, where 



/ = Jtr|VQ| 2 + ^trQA B + i;_(trQA E ) 2 + ^trQA 2 , 



A# = fj 3 (g) 1, X z = 1 ® ct 3 , u_ 



E : 
~i 

NB 



V- (12) 



Here, J is given by Eq. pTj) up to numerical prefactor of 
order unity, and we recognize tr QXe, and trQA 2 to be 
the degrees of layer and spin polarization respectively. 
Linear B dependence of the 'capacitance energy' term 
arises because the density of electrons in the zeroth Lan- 
dau level is proportional to B. The Zeeman term is rel- 
atively small, Ez/v^ sa (m/m e )/€ <C 1, where m c is the 



free electron mass, and can thus be neglected. The ca- 
pacitance term, meanwhile is parametrically smaller than 
stiffness, u_/J w d/ao <C 1, which ensures validity of the 
sigma model for not too strong electric fields. 

The sigma model, Eq. lfT^)) . can be used to assess sta- 
bility of the QHFM state in the presence of electric field. 
The Hamiltonian, Eq.(j4|), at E ^ has symmetry group 
SU(2) S pi n ® U(l) va uey, broken by the QHFM state for any 
choice of Q other than Q = 1 (g> ffy. In the absence of 
capacitance and Zeeman energies, Eq. (fT2"]) predicts the 
state Q — 1 <£> 773 for any nonzero E. However, for small 
electric fields, such a state is disfavored by capacitor en- 
ergy. Reinstating it in Eq. (fT2|) . we see that the symmetry 
preserving state is only realised for e > e C 2 = NV-B /$q. 
This critical field marks a second order phase transition 
from QHFM state to layer polarized state [diagonal line 
in FigQJa)]. We note that since e C 2 is less than J by a 
factor of about d/ao, the phase transition occurs within 
the realm of validity of the sigma model. The linear de- 
pendence s C 2 oc B is in agreement with experiment Q- 

In summary, we have argued that at small E and B 
BLG hosts an SU(4) flavor singlet AHI state which breaks 
time reversal symmetry, and exhibits QHE even at zero 
magnetic field. We have shown that the application of 
an electric field restores time reversal symmetry, while 
application of a magnetic field induces SU(4) symmetry 
breaking transition to QHFM state. We have calculated 
the properties of the QHFM state, taking full account 
of screening, and have predicted an electric field induced 
phase transition to a time reversal invariant state. 
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